MEASURING AGGLOMERATION OF AGGLOMERATED PARTICLES PICTURES 
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Abstract. In this article, we introduce a novel geometrical index 5agg, which is associated with the Euler number and is 
obtained by an image processing procedure for a given digital picture of aggregated particles such that 6agg exhibits the 
degree of the agglomerations of the particles. In the previous work (Matsutani, Shimosako, Wang, Appl. Math. Modeling 37 
(2013), 4007-4022), we proposed an algorithm to construct a picture of agglomerated particles as a Monte-Carlo simulation 
whose agglomeration degree is controlled by 7agg G (0, 1). By applying the image processing procedure to the pictures of the 
agglomeration particles constructed following the algorithm, we show that (5agg statistically reproduces the agglomeration 
parameter 7agg. Keywords: agglomeration, digital image processing procedure, Euler number 
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1. Introduction 

Nano-composite materials have a promising future from 
industrial viewpoints, since in the materials, geometrical 
properties in micro-scale play crucial roles and generate 
novel and various macro-material properties. By control- 
ling the geometrical properties or shapes, we can design the 
macro-material properties drastically. Following Kelvin's 
philosophy of science [s]^, it is quite important to evaluate 
such geometrical properties or shapes if one needs to control 
them using the scientific knowledge. 

On the other hand, the smaller the particle is, the larger 
the effect of the surface energy is. It means that small parti- 
cles are apt to aggregate or agglomerate in general because 
the agglomeration and the aggregation of the particles de- 
crease the total surface energy and contribute to the stabil- 
ity of the system. When we handle materials consisting of 
nano-particles, the agglomeration and the aggregation are 
ones of the most important shapes since they sometimes 
have an effect on the generations of the macro-materials 
properties. In this article, we focus on them. It is said 
that the aggregation is due to chemical effects whereas the 
agglomeration comes from physical effects. Since in a com- 
putational model, there is no difference between them, we 
call both agglomeration in this article, though in spatial 
point analysis [6 , the aggregation is chosen in general. 



^ Kelvin wrote his philosophy of science, "In physical science the 
first essential step in the direction of learning any subject is to find 
principles of numerical reckoning and practicable methods for measur- 
ing some quality connected with it. I often say that when you can 
measure what you are speaking about, and express it in numbers, you 
know something about it; but when you cannot measure it, when you 
cannot express it in numbers, your knowledge is of a meagre and un- 
satisfactory kind; it may be the beginning of knowledge, but you have 
scarcely in your thoughts advanced to the state of Science, whatever 
the matter may be." ^ 



In the article [10 , in order to find the agglomeration ef- 
fect in the electric conductivity of the nano-composite ma- 
terial, we study the electric conductivity in an agglomerated 
continuum percolation model and show that the agglomera- 
tion of particles affects the macro-material properties. The 
purpose of this article is to evaluate the agglomeration in 
the binary digital images of agglomerated particles, e.g., of 
electron-microscopes. 

For the same purpose, so many evaluation methods and 
definitions of the agglomeration are proposed to evaluate the 
agglomeration. In spatial point analysis, the distribution of 
the nearest distance particles, Clark-Evans index and so on 
are considered [2^, ^6 . Further Miles considered the problem 
in [11] and showed the two-dimensional overlapping ratio 
of the random configurations in three dimensional space. 
These investigations on the agglomeration have been done 
in the framework of the statistical analysis for a point pat- 
tern 1Z = {pi e M?} which are given as statistical configu- 
rations of (finite) points. In the analysis, IZr := UpG7^ ^^^p 
is investigated, which is a configuration of disks whose cen- 
ters are 1Z, where Us,p := {q e R'^ \ \q - p\ < e }. The 
Euler number, the area and the perimeter of IZr for several 
point processes 7^'s are computed as morphological indeces 
or the Minkowski characterization [6 . When TZ is given by 
the point process of the Poisson type, Stoyan, Kendall and 
Mecke studied their behaviors based on the study of Miles 
pT] and found that 

(1) 

e{x) = (1 - x)e~^, a{x) -(1 - e~^), l{x) e~^. 



where e(x), a(x), and ^{x) are the normalized versions of 
the Euler number, the area and the perimeter of T^r, and x 
is a normalized radius r [HI [8] . Mecke and Stoyan studied 
the difference among point patterns given by different pro- 
cesses in terms of these behaviors 1 8 . Further Tscheschel 
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and Stoyan also studied the statistical reconstruction of ran- 
dom point patterns |T5] . 

However in the nano-composite materials consisting of 
nano-particles, the particles themselves sometimes have com- 
plicated shapes, such as ellipsoids and rods, as we investi- 
gated in [9 . In other words, the configuration is not given 
by a point pattern with radius r in general and further its 
parts are overlapped like (d) in Figures [2] and [3l Hence it is 
basically an ill-posed problem to define the center points of 
actual agglomerated particles in a given picture, e.g., of an 
electron-microscopes. 

Thus it is natural to consider geometrical properties of 
the binary picture as a general geometrical object embedded 
in R^. 

Recently MacPherson and Schweinhart proposed a novel 
method which evaluate the complicatedness of the compli- 
cated geometric objects embedded in a plane M? in terms of 
the persistent homology [7 . The persistent homology gives 
the homological quantities of the filtered persistent modules 
with real parameter [3l [16] . It could be regarded as a gen- 
eralization of homotopical approach in traditional algebraic 
topology [1 , though the deformation does not preserve ho- 
motopical properties. For a geometrical object M C M^, we 
consider a family of objects with a real parameter t G [0, 1], 
i.e., {Mt \te [0,1]}. 

By considering union of the e-neighborhood of each point 
in A^, Ms '= [jp^M ^^'P' induced from the standard Eu- 
cHdean topology, MacPherson and Schweinhart evaluated 
the complexity of the geometrical objects. The persistent 
homology shows the distributions of topological changes 
generated by the filtered persistent modules (vector spaces) 
induced from A4t' C A4t for f < t. 

As in [7 , to investigate the effect from the standard 
topology of Euclidean space and to evaluate the complexity, 
we use the one-parameter family of a deformed geometrical 
object, and propose a digital image processing procedure 
which characterizes the shapes in pictures of the electric 
microscope in this article. (In Section O we give the list of 
assumed geometrical features of the pictures which we deal 
with.) For an appropriate geometrical object A4 in with 
a characteristic lengths ii and £2 , we also handle the family 
of geometrical objects {Mt = [jpeM ^t,p \t G [h^h]}- We 
define the cumulus of the absolute differential Euler number 



(CADE) by, 
(2) 



-I" 



dx{Mt) 



dt 



dt 



where x{^) is the Euler number of X. S{M] i2^£i) evalu- 
ates how many topology changes occur for the deformation 

[eui2]. 

As we are concerned with the image processing procedure 
for images of the electron-microscopes, we will customize 



S{M;i2^h) as S{M;i2^h) for any binary pictures as an 
image processing procedure, which is shown in Section [3] 
more precisely. Further in the nano-materials, there are 
several scales and one of them is the size of the particles and 
the resolution of the digital picture is given by the pixel size. 
We fix £1 and £2 by the pixel size and the (average) radius 
p of the particle respectively to evaluate the agglomeration 
and propose an agglomeration index, 

(3) (^agg(A^) := j^i^piM) -S{M;p,a)), 

where p{M) is the volume fraction of M in the region W 
{M (ZW (Z M^), 8p is the average of a "standard pattern 
of volume fraction p" as mentioned in Section [3l and a is a 
normalized factor 1.2. 

In order to estimate our agglomeration parameter (^agg, 
we performed Monte-Carlo simulations for the binary ag- 
glomeration configurations of particles whose degree of the 
agglomeration is parameterized by 7agg, since in the arti- 
cle [10], we proposed a statistical model which numerically 
generates the agglomeration of particles controlled by the 
parameter 7agg in order to investigate the properties of the 
agglomerated continuum percolation models. In Section [21 
we review the algorithm following the article ^lOj. For a 
given parameter 7agg G [0, 1], we can statistically construct 
the infinitely many configurations with the same level of the 
agglomerations. As the continuum percolation model is the 
same as a germ-grain model in the study for the point pro- 
cess [6 , there are several other algorithms to construct ag- 
gregational germ-grain models, such as Neyman-Scott pro- 
cesses [6 , though they are different from ours; ours is for the 
actual pictures of electron-microscopes of the agglomerated 
nano-composite materials with the radius p as mentioned in 
[10 . We apply the index (5agg to evaluate the agglomeration 
of the agglomerated configuration which is generated by the 
forward method in [TO]. Then the relevancy between ^agg 
and 7agg is shown in Section IH i.e., in Figure [71 and Table 
[3l We also mention the relation among our (^agg and the 
well-established Clark- Evans index in Section IH 

2. Agglomerate configuration 

In order to explain what is the agglomeration that we are 
concerned with, we show the agglomeration configurations 
in computer science, which we handled in [10]. In the article 
p[Q] , we proposed a construction of the agglomerated contin- 
uum percolation models which apparently recover geometric 
properties of real nano-particles, though there are several 
other agglomerated percolation models such as Neyman- 
Scott processes [6] |14| |15]. Since our method has a only 
single parameter 7agg besides a typical length p whereas 
others are given as a point process with several parame- 
ters, we believe that ours is a natural model for the actual 
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pictures of electron-microscopes of the agglomerated nano- 
composite materials with the radius p. As shown in Figures 
[2] and [31 we have the agglomerated configuration of particles 
depending upon a agglomeration parameter 7agg G [0, 1]. In 
this section, we show the geometrical setting of agglomer- 
ated continuum percolation model in [10 , which is modeled 
by the agglomerated clusters in nature. 

We set particles parameterized by their center positions 
{x^y) into a box-region W := [0,L] x [0,L] at random and 
get a configuration A^n as a model of continuum perco- 
lation. The particle corresponds to a disk with the same 
radius p, ^a^,,^, := {{x,y) G W | \{x,y) - {xi,yi)\ < p}. The 
configuration Mn is given by Mn '= [j^=i ^Xi^yi- 




Set a disk at a random position (x,v) 




Figure 1 . The fiowchart of the algorithm 
which constructs the agglomerated config- 
urations. 

The fiowchart in Figure [T] illustrates the algorithm. As 
the initial state, the configuration Mo has no particle. As 
the first step, for a uniform random position (x, y) G W, we 
set a particle Bx^y whose center is (x, y) and the radius is 
p, i.e.. Ml := Bx,y. 

For the (n + l)-th step, we take a position (x, y) at uni- 
form random in W, and another random parameter 7 at 
uniform random in [0, 1]. If 7 is greater than 7agg, we set 
Mn-\-i •= Mn[j Bx^y. We now allow the particles to overlap 
each other. 

For the case 7 < 7agg, we first check whether the disk 
Bx,y is connected with the previous configuration Mn or 
not. For the case Mn H ^x,y 7^ 0, we employ the position 
and set Mn-\-i '= Mn[j Bx^y. Otherwise or Mnr\^x,y = 
0, we abandon the position and go on to take another uni- 
formly random position (x, y) in W until we find the position 
which supplies a connected particle Bx^y with Mn- 

In other words, for the case 7 < 7agg, the added parti- 
cle must be connected with the previous configuration Mn- 



Thus, 7agg stands for the agglomeration of the particle sys- 
tem. 

By monitoring the total volume fraction which is a func- 
tion of Mn and is denoted by p{Mn)^ we go on to put the 
particles as long as p{Mn) < p for a given volume frac- 
tion p. We find the step n{p) such that p{Mn{p)-i) < P 
and p{Mn{p)) > P- Since we assume that the difference 
between p{Mn{p)-i) and p{Mn{p)) is sufficiently small, we 
regard p{Mn{p)) as the volume fraction p itself hereafter 
under this accuracy. 

Since in the Monte-Carlo method, we use the pseudo- 
randomness to simulate the random configuration Mn{p) 
for given p and 7agg, the configuration Mn{p) depends upon 
the seed is of the pseudo-randomness which we choose. We 
let it be denoted by M^^^^^p^ig or its statistical quantity by 




(a) (b) 




(c) (d) 

Figure 2. The agglomerated configura- 
tions ofp = 0.2: These (a), (b), (c), and (d) 
show the configurations with the agglomer- 
ation parameter 7agg = 0.0, 0.3, 0.6 and 0.9 
respectively. 

For sufficiently large L and V {p <^ L' < L) and for a 
window := [x, x + L'] x [y, y L^] C W, the volume 

fraction p'{x^y) in W^'^ is proportional to the number of 
particles N'{x^y) in W(^^), i.e., there is a constant num- 
ber A such that p'{x^y) = \N'{x^y). Further M^^^^^p^i^ is 
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Figure 3. The agglomerated configura- 
tions ofp = 0.4: Tiiese (a), (b), (c), and (d) 
siiow ttie configurations witii ttie agglomer- 
ation parameter 7agg = 0.0, 0.3, 0.6 and 0.9 
respectively. 



isotropic and independently scattered and thus the asymp- 
totic behavior of A^7agg,p,is is a kind of Poisson process [6l 
p. 66] though we have the condition L' ^ p. 

Further though our algorithm for 7agg > is not Markov 
process because each step depends on the previous configu- 
ration, we should note that it preserves a hierarchical struc- 
ture for fixed is and p > p' ^ 

As we will show the assumed geometric properties of pic- 
tures in Section [3l the pictures which we will deal with are 
illustrated in Figures [2] and [3l which are obtained follow- 
ing our algorithm. Our more concrete aim of this article 
is to recover the 7agg for a given configuration M-^^^^^p in 
a statistical meaning. More precisely, our study is to find 
statistical monotone functions of 7agg and to show that one 
of them is ^agg in (|3|). 

3. Evaluation of Agglomeration for a 
configuration mp^^^^^ c w 

As we are concerned with the evaluation method as a 
digital image processing procedure [12^, in this section, we 



illustrate our algorithm for a picture which only has binary 
values. It is natural that we assume the configuration M. 
(implicitly A^p,7agg,i5 and a picture of nano-composite mate- 
rial in an electron- microscope) has the following structures: 

(1) L is sufficiently larger than p so that the particles 
of M. are a representative of sufficiently random- 
ized configurations; we could assume the Euclidean 
invariance (translation, rotation and inversion) sta- 
tistically; after averaging them, the physical and ge- 
ometrical quantities are invariant for any Euclidean 
action E(2) up to the statistical deviation. If the 
deviation is not small, we could consider the series 
of {A^p,7agg,is \ '^s}' (It means that for the case of 
the pictures of the electron-microscopes, we could 
assume that the researchers prepare the series of 
pictures of a material or materials which are pro- 
duced in the same conditions.) 

(2) It is assumed that the volume fraction is less than 
the percolation threshold of two dimensional con- 
tinuum percolation models. (For the case of nano- 
composite material which is based upon the perco- 
lation theory, the volume fraction around the per- 
colation threshold 0.2 ~ 0.3 in three dimensional 
percolation models is concerned, which is far less 
than the percolation threshold of two dimensional 
case 0.5 0.7.) 

(3) There are three sizes of the system or the picture 
•^P,7agg (^nd digital image of nano-composite 
material of an electron-microscope); 

(a) the (average) size of particles, which is given 
by P, 

(b) the analyzed size of the system, which is, now, 
given by L as mentioned above, and 

(c) the pixel size a, which is also controlled so that 
we can discriminate the particles in concerned 
resolution. 

Under these assumptions, we consider geometry of A^. It 
is known that the e-neighborhood, M.^ = U^gM ^^^p H 
can be realized by the so-called level set method in com- 
puter science [13 . Let d : W ^ M be the signed distance 
from the boundary dM so that the outer side is assigned 
to the positive distance and the inner side is to the nega- 
tive one, and then the geometrical object in the level set 
method can be regarded as Ct = d~^{t). A4t of {t > 0) is 
equal to d~^{[0,t])[j M and Ct = dMt- For t < case, 
= ^{[jpedM ^t,p) \ C-t- Hence by means of the level 
set method, we can compute the more precise geometrical 
properties beyond the pixel size resolution even on the im- 
age defined over a subset of I? . 
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However in the digital image processing procedure, we 
investigate the geometrical object up to the pixel size res- 
olution in general. Further we must pay our attentions on 
the computational cost if we apply our method to real prob- 
lems in industry, though level set function method requires 
higher computational cost than a simple digital image pro- 
cessing procedure. Hence in this article, we use the thicken- 
ing scheme in the image processing procedure [12] instead 
of the level set function. Though the ordinary thickening 
scheme has anisotropic behavior, it does not have a serious 
effect on the result because the configuration itself is iso- 
topic or rotational invariant. We use the modified thicken- 
ing scheme, which improves the anisotropic behavior shown 
in Section m Let A^*^*^ be the i-th thickening of M in W. 
We modify the CADE ([2]) as an image processing procedure 

by 

i=ni + l 

In the persistent homology, the Betti number is handled 
in general. Since the computational cost to the evaluation 
of the Euler number is not so high and the Euler number 
could be compared with the results in [6, 14 , we consider 
the behavior of the Euler numbers of A4t in this article. 
More precisely though there is no guaranteed that xi-^t) 
is equal to x{M^'~^^) for t G [a{i - 1/2), a(i + 1/2)), we 
handle x(-^^*~^^)5 mentioned above, in digital analysis, 
we should basically neglect finer geometrical difference than 
the pixel size resolution and we follow the principle. In the 
complicated system, we believe that it is quite important 
how many topology changes occur for the z-step, and the 
difference of the Euler number can represent the behavior. 

Further the agglomeration can be discriminated whether 
the particles are connected or not. From ([T]) , if L/p is suffi- 
ciently large, even for 7agg = and small S{A4o^p^na^O) 
does not vanishes n > 0, in general, due to the randomness 
of the configurations. Further the behavior S{A4;nia^n2a) 
of ni,n2 G [0,p/a) is quite important since the agglom- 
eration suppresses the topology change in the interval as 
illustrated in Figures H] and [51 

Due to the randomness of the configurations and the ag- 
glomeration, it is not so important whether the Euler num- 
bers increase or decrease, but the topological change for the 
deformation is quite important. We define the agglomera- 
tion parameter (^i^^'^^^ in dU more precisely 

(4) ^i-'-)(A^) = 

PiM) 

where p{Ai) is the volume fraction of A^, ^^^^'^^^ is the av- 
erage of the standard patterns of volume fraction p, and a 
is a normalized factor 1.2, which is chosen as a result of the 



comparison with 7agg (see Table [3]). The standard pattern 
means the pattern of 7agg = with the same radius in the 

same window W. Then ^agg'^^^(A^7agg,p) characterizes how 
many topological changes occurs in the interval (nia, n2, a) 
for the deformations for each particle in M^^^^^p by normal- 
ized by 4"""'^. 



4. Numerical Computation and Results 

Let us show the relevance between Sp and jp by the 
Monte- Carlo simulations following the algorithm mentioned 
in Section [21 Using the algorithm in Section [21 we have ten 
pictures of agglomerated particles for each 7agg = 0, 0.3, 0.6 
and 0.9, and for each p = 0.1,0.2,0.3 and 0.4 by letting 
L = 2400 and p = 10 as in Figures [21 and [3l 

On the thickening to compute the CADE, we use two 
types thickening process, 

□ 

type I: □ ^ □ □ □, 

□ 

□ □ □ 
type II: □ ^ □ □ □, 

□ □ □ 



such that we generate an octagon asymptotically and ap- 
proximates the area of the disks; In other words, on the 



Table 1. The pattern of thickening 



steps 


type 


radius 


area 


n.of pixels 


1 


II 


0.5 


0.785 


1 


2 


I 


1.5 


7.065 


9 


3 


I 


2.5 


19.625 


21 


4 


I 


3.5 


38.465 


37 


5 


II 


4.5 


63.585 


69 


6 


I 


5.5 


94.985 


97 


7 


I 


6.5 


132.665 


129 


8 


II 


7.5 


176.625 


185 


9 


I 


8.5 


226.865 


229 


10 


I 


9.5 


283.385 


277 



thickening process, we use the deformation in digital pro- 
cess procedure for each point which is given in Table [H 
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Further for each point, we consider the thickening: 



We computed the ten pictures for eachp = 0.1, 0.2, 0.3, 0.4 
and 7agg = 0.0, 0.3, 0.6, 0.9 with ten random seeds. Figure S] 
shows the CADE, f(A^p,7agg; '^cl^ 0)^ foi" the n-th thickening 
step for each p. 
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Figure 4. The CADE S{Mp,j^^^;na,0) 
of the n-th thickening step for each 
7agg; Those of the volume fraction p = 
0.1,0.2,0.3 and 0.4 are ihustrated in (a), 
(b), (c), and (d) respectively. 

On the other hand, though it is difficult to identify the 
center points of the particles for given pictures, especially 
of the agglomerated case as shown in images (b) and (c) of 
Figures [2] and [31 we know the data of the center points of the 
particles. Thus we can use the techniques of the statistical 
analysis for the spatial point patterns. Figure [5] displays the 
global distribution of the Euler numbers of different radius 
of a seed by using the software provided in [6, p. 2040. 



Since our radius is 10 and our agglomeration algorithm is 
characterized by the radius, the behavior of the distribution 
in Figure [5] strongly depends on the regions p > 10 and 
p < 10. Figure m correspond to (10,20] region and thus it 
implies that our improved thickening algorithm works well 
except the first thickening step of the 7agg = 0.9. 

Since the agglomeration in our algorithm means that the 
number of agglomerated particles is larger than the uniform 
randomness 7agg = 0. The variation of the Euler number 
is related to the deformation in which disjoint clusters con- 
nect due to the thickening. Agglomeration means that the 
number of the disjoint clusters is less than that of uniform 
randomness. The variation of the Euler number for the in- 
creasing of the radius p > 10 in Figure [5] is suppressed for 
large 7agg- Hence the dependence in Figure |4] is very natural 
except the first thickening step of the 7agg = 0.9. 

Further in the image processing procedure, we must pay 
attention to the digit alized errors. We should recognize that 
the first step contains some digitalized errors because the 
behavior in Figure |4] is contradict with that in Figure [5l the 
behavior of the curves of 7agg = 0.9 in Figure [5] are very 
mild over (10,20) whereas the first steps of 7agg = 0.9 in 
Figure m rapidly increase. 

Hence we are concerned with ^(A^p,7agg; p, «). 
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Figure 5. The Euler number vs the ra- 
dius p as the point pattern of M^^^^^p^ig'. 
(a), (b), (c), and (d) illustrate the Eu- 
ler numbers of the volume fraction p = 
0.1,0.2,0.3 and 0.4 respectively. 

Tableland Figure [6] show the dependence of the CADE, 
£^(A1p,7g^gg; p, a), on the agglomeration parameter 7agg for 
each p. They exhibit the negative correlations. The agglom- 
eration means that the number of agglomerated particles is 



Table 2. CADE vs 7agg: 
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Figure 6. The CADE, S{Mp^^^^^; p,a), 
vs 7agg: (a), (b), (c), and (d) display 
the states of the volume fraction p = 
0.1,0.2,0.3 and 0.4 respectively. 
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Figure 7. The agglomeration index 
^agg VS 7agg: (a), (b), (c), and (d) dis- 
play the states of the volume fraction p 
= 0.1,0.2,0.3 and 0.4 respectively. 

larger than the uniform randomness 7agg = as mentioned 
above. Since the agglomeration prevents the topological 
changes on the thickening, we have the negative correlation 
in Figure [6l 

We, now, define the the average £p of a "standard pat- 
tern of volume fraction p" by the average of the CADE, 



p 


0.1 


0.2 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





805.7 


879 


750 


2131.9 


2221 


2043 


0.3 


582.9 


657 


514 


1501.2 


1675 


1388 


0.6 


317.6 


340 


269 


821.9 


886 


766 


0.9 


141.9 


204 


78 


244.0 


364 


160 


P 


0.3 


0.4 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





2878.5 


2954 


2806 


2705.3 


3072 


2555 


0.3 


2035.2 


2181 


1910 


1966.3 


2160 


1798 


0.6 


1121.9 


1262 


1003 


1066.6 


1188 


905 


0.9 


408.8 


507 


248 


582.1 


702 


464 



^(^p,o,is; P^ci)^ of the uniform random configuration, i.e.. 



Tagg — case, 



1 ^° 

25 = 1 

Thus we denote the agglomeration index ^igg^'"*^^ by ^agg as 
in (|3j). Further in order that Jagg corresponds to 7agg, we 
chose a = 1.2. 

Figure [3 and Table [3] show the relation between Jagg and 
7agg; both show that 7agg is correlated to Jagg and approxi- 
mately recovers Jagg up to the statistical fiuctuation. 

Table 3. Jagg vs 7agg: 



p 


0.1 


0.2 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





0.000 


0.083 


-0.109 


0.000 


0.050 


-0.050 


0.3 


0.332 


0.434 


0.221 


0.355 


0.419 


0.257 


0.6 


0.727 


0.799 


0.694 


0.737 


0.769 


0.701 


0.9 


0.989 


1.084 


0.896 


1.063 


1.110 


0.995 


P 


0.3 


0.4 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





0.000 


0.030 


-0.031 


0.000 


0.067 


-0.163 


0.3 


0.352 


0.404 


0.291 


0.328 


0.402 


0.242 


0.6 


0.732 


0.782 


0.674 


0.727 


0.799 


0.673 


0.9 


1.030 


1.097 


0.989 


0.942 


0.994 


0.889 



In the statistical analysis of the spatial point patterns, 
the Clark- Evans index is a well-established index which rep- 
resents the agglomeration degree of a given point pattern, 
though in general, it is very difficult to identify the center 
points of the particles for a given picture, such as images 
(b) and (c) of Figures [2] and [3l the problem is sometimes 
ill-posed for the cases. Since we know the data of the cen- 
ter points of the particles of every A4^^^ 



we illustrated 

the Clark-Evans index in Table H] and Figure [HI which show 
that the Clark-Evans index represents our agglomeration 



8 



parameter 7agg well. The correlation between the Clark- 
Evans index and S^gg is displayed in Figure [H It shows a 
good negative-correlation for each volume fraction p. 
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Figure 8. The Clark-Evans index 
vs 7agg: (a), (b), (c), and (d) display 
the states of the volume fraction p = 
0.1,0.2,0.3 and 0.4 respectively. 



Table 4. Clark Evans index vs 7agg- 



p 


0.1 


0.2 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





1.018 


1.030 


0.999 


1.011 


1.022 


0.995 


0.3 


0.755 


0.776 


0.733 


0.837 


0.850 


0.828 


0.6 


0.574 


0.599 


0.561 


0.703 


0.709 


0.693 


0.9 


0.412 


0.419 


0.402 


0.556 


0.569 


0.547 


P 


0.3 


0.4 


7agg 


Ave 


Max 


Min 


Ave 


Max 


Min 





1.006 


1.010 


1.002 


1.003 


1.008 


1.000 


0.3 


0.890 


0.897 


0.881 


0.927 


0.933 


0.922 


0.6 


0.788 


0.795 


0.783 


0.850 


0.856 


0.845 


0.9 


0.661 


0.665 


0.655 


0.741 


0.746 


0.735 



5. Summary 

In this article, we introduced the novel geometrical index 
<^agg, which is associated with the Euler number and is ob- 
tained by an image processing procedure for a given digital 
picture of aggregated particles such that ^agg represents the 
degree of the agglomerations of the particles. Following the 



Figure 9. The Clark-Evans index and Jagg- 

algorithm in [10 , we constructed digital pictures of aggre- 
gated particles controlled by the agglomeration parameter 
7agg G (0, 1) as a Monte-Carlo simulation. By applying 
the image processing procedure to the pictures, we showed 
that Jagg statistically reproduces 7agg. Since we have the 
data of the center points of the particles, we also computed 
the well-established Clark-Evans index and showed that it 
also represents 7agg well. However though the methods in 
the point process analysis including the Clark-Evans index 
require the data of the configuration of the points, the de- 
termination of the center points of the particles for a given 
picture is basically an ill-posed problem. Hence our method 
has an advantage because we don't need to find the center 
points in the computation of Jagg- In other words, our pur- 
pose that we recover 7agg for a given picture by means of 
the digital image processing procedure is accomplished by 
considering the deformation of the geometrical object. It 
implies that we can measure the agglomeration in a given 
picture of agglomerated particles. 

In this article, we have investigated pictures whose vol- 
ume fraction p is less than 0.5, because it is difficult to deal 
with pictures with large volume fraction. It is expected to 
find further natural index to discriminate the agglomera- 
tion with the large volume fraction, e.g., in terms of the 
persistent homology [7]- 
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